Figure Captions 



Fig.l. Region of absolute convergence of the solutions J3 and J4. 

Fig. 2. Region of absolute convergence of the H2 hypergeometric function. 

Fig. 3. Region of absolute convergence of J2[l], and Je- The curve that 
separates J2 of the others is a branch cut, see eq.(46). 

Fig. 4. Region of absolute convergence of the Jg. The solution Jj is symmetric 
to it in s ^ t. They are finite and hold in the relativistic regime of forward 
scattering in the t and s-channel respectively. 

Fig. 5. The general 2-particle — >■ 2-particle Feynman graph that has leading 
singularity (ies) on the physical sheet [17]. 
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Abstract 



In this sequel calculation of the one-loop Feynman integral pertaining to a massive 
box diagram contributing to the photon-photon scattering amplitude in quantum 
electrodynamics, we present the six solutions as yet unknown in the literature. These 
! six new solutions arise quite naturally in the context of negative dimensional 

integration approach, revealing a promising technique to handle Feynman integrals. 
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1 Introduction. 



In our previous paper[l] we calculated the Feynman integral pertaining to the 
photon-photon scattering amplitude in Quantum Electrodynamics (QED) us- 
ing a technique we call negative dimensional integration method (NDIM)[2-5]. 
One of the outstanding features of NDIM is that the complexities of perform- 
ing D-dimensional momentum integrals are transferred to the easier task of 
solving systems of linear algebraic equations. In [1] NDIM has allowed us to 
recover very easily the two known hypergeometric series representations for 
the pertinent Feynman integral. Another outstanding feature of NDIM — and 
this is in our opinion its greatest potential — is that it gives simultaneously 
six new results for the integral in question in a very straightforward manner. 
The aim of the present paper is to consider these new results. Each of them 
is valid in certain regions of external momenta and are related to the others 
by analytic continuation, either directly or indirectly. 

Compared to the traditional methods of working out Feynman integrals, NDIM 
is by far simpler and more straightforward. The difficulties of calculating para- 
metric integrals — if one is handling a Feynman integral in the standard way[6] 
— or solving contour integrals — if one is using the Mellin-Barncs' integral 
representation for massive propagators [7-9] — are with finesse embedded in 
an easier problem of solving systems of linear algebraic equations. 

The Mellin-Barnes' integral technique can provide us with only two types 
of results: cither the resulting function will depend on, say, p^/m^, or its 
inverse, m^/p^, depending on whether one chooses to close the contour to 
the left or to the right in the complex plane. The resulting functions arc in 
general hypergeometric ones. However, if we have more than one external 
momentum (as in our present case) and/or internal mass, there are several 
other combinations of dimensionless variables that define distinct regions of 
external momenta. It is then clear that the Mellin-Barncs' technique will be 
limited to give simultaneously only two among those several possibilities. Not 
so with NDIM, where many of the possibilities are accounted for at the same 
time. Of course, in both cases, one can always construct other power series 
representations by means of suitable analytic continuations as long as the 
formula for such extensions be known. NDIM has then clear advantages over 
other techniques in this sense. 

Hypergeometric functions of one and two variables have lots of well-known 
analytic continuation formulas but as the number of variables increases — as 
far as we know — the fewer the known rclations[10] are. On the other hand, 
since NDIM provides us with very many simultaneous results, which in prin- 
ciple must be connected by analytic continuation, we come to the realization 
that it is not only a very powerful technique to work out Feynman integrals 
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but an elegant approach to check on analytic continuation properties of the re- 
sulting functions as well. Consider our present case: Altogether we have eigth 
distinct solutions for the Feynman integral for the photon-photon scattering, 
two of which have already been considered in our previous paper and six new 
ones, which are connected with each other by suitable analytic continuation 
formulas. 

One could rightfully ask: Why do we want so many distinct results at the 
same time? We give some good arguments for this. Firstly, if we have only 
oneQ or at most two[^ results in distinct regions of the external momenta, 
all the other regions must be worked out all the way through the analytic 
continuation formulas, which is not always an easy task to perform and cer- 
tainly is very much time-consuming. Secondly, the important special case of 
forward scattering in the relativistic regime cannot be dealt with if one has 
only the two known hypergeometric series representations for the Feynman 
integral relative to the photon-photon scattering. These series are unsuitable 
for handling this special case because of the very nature of their variables. 
The same reasoning applies to the backward scattering. Thirdly, our results 
are expressed in a compact form that can be transformed — by an appropriate 
integral representation — into the more cumbersome standard form in terms 
of dilogarithms, if one wants to do so. Fourthly, we can identify the branch 
points and singularities of Feynman integrals directly from their hypergeomet- 
ric series representations. Fifthly and lastly, since any two distinct solutions 
are related by analytic continuation, NDIM is an elegant and economical way 
of obtaining analytic continuation formulas among hypergeometric series (see 
Section 3). 

The outline for this paper is as follows: In Section 2 we write down the new 
results for the Feynman integral in D-dimensions and arbitrary powers of 
propagators much the same as in dimensional regularization[ll]. In all of the 
resulting expressions we have singularities (poles) either when we want to take 
the physical limit of D = 4 or in the limits of unity powers for propagators, 
so that Section 3 is devoted to analyse the special cases we are interested 
in by introducing suitable regularization parameters. Four of the solutions re- 
main divergent even after the regularization procedure and deserve our special 
analysis and discussion. We also examine the convergence regions for the new 
results and draw diagrams for all of them. Finally in the last section. Section 
4, we make our concluding remarks and present some further challenges and 
applications for NDIM. 



^ via parametric integration 
^ via Mellin-Barnes' integration 
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2 New Results from NDIM. 



In our previous paper[l] we studied the integral, 



J{i,j,k,l;m) = j d^q [q^ - 



m 



{q — pY — {q — k\Y — w? 



X 



{q - k^f - 



(1) 



whose counterpart in positive D is 



dPq 



{q^ — rri^y [{q — p)^ — m?Y [{q — ki)^ — rm?] 
1 



X- 



[(? - ^2)^ - rn?] 



which is relevant to the photon-photon scattering in QED. In particular, we 
are interested in the case J(— 1, —1, —1, —1; m) = K{1, 1, 1, 1; m). 

Hereafter, unless otherwise noted, we closely follow the notation of [1]. So, let 
a — i -\- j -\- k -\- I -\- \D , where D is the space-time dimension and let the 
Pochhammer symbol be 



[a\k) = {a)k 



T{a + k) 



Here we remind ourselves that the results of (1) come from the solutions of a 
system of linear algebraic equations and that for this particular integral the 
system is such that has twenty-one solutions altogether, six of which are trivial 
ones and fifteen of which give double series. In [1] we listed five of them which 
combined appropriately among themselves led to the two well-known hyper- 
geometric series representations [8] of the Feynman integral. Below we give the 
remaining ten solutions for the system, which combined appropriately among 
themselves, yield the six new solutions for the pertinent Feynman integral. In 
the following we use s and t for the usual Mandelstam variables. 



h = /e Si yae, ag, Pe, /3g, Oq; 73, 6lg 

It = iQ^i ^ k, j ^ I \ s ^ t) 

h = /s as, /?8, 78, ^8, (ps; Ps, (p's 



-t 4m^' 



-t -4m^' 



t 



(2) 

(3) 
(4) 
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Ig^Is{k ^ I) 

ho = fio S2 ^ttio, Pio, 7io, ^10, 0io; Pio, 0'io 
hi^hoii ^ 3) 

lu = /12 S2 (^ai2, P12, 7l2, (^12, 012; P12, 0'l2 
Il3 = Il2{k ^ I) 

lu — fu S2 ^Q;i4, Pu, 7i4, 5i4, 014; pu, 0'i4 
i"l5 = -^14(^ ^ i) , 



—s — 4m^' 
4m^ -t \ 



s ' 4m2^ 
4m^ —s ' 



t '4m2 



where we have defined the following two functions 



«Si(Q;,Q;,/i/,/^,6';7,e' |2;i,2;2)= 2^ — j-j- — 



X 



and 



S2[a,(i,-f,d,(p;p,(p\zi,Z2)^ 2_ — j-r 



(pI/^) 



X 



(0k -/^) 



with 



/6 = 

/lO = 
/l2 = 
/l4 = 



^ V4 



(-fc|-^-J-|D)(-' 


-/.■ + /) (-i 


a-/) 


{—i — / + 0" p 





-,■ (-^1 - ^ - J - + - - fc) 



(-i - A; + cr|i + /c + ^D) 
-fc|0(a+|L'|/-2a-|D) 

-%\2){o^\D\3-2a-\D) 



with the following parameters for the function Si. 

^' 1 1^ In 

"6 = 2 ~ 2*^" 4-^' 

e, = -a, 

76 = 1 + i - 
and the parameters for the function S2'- 





a8^j + k + \D, 


QJiO = 






Ps = -I, 


/3io = 


-J, 




78 = 1 - ifT - \D, 


7io = 






= 1 - - i^, 


(^10 = 






(f)s^ -i-j -k - \D, 


010 = 


-i-k-l-\D, 




Ps^l + k-l, 


PlO = 


i + i-3, 




(p'^ = l-i-k-\D, 


0'io = 


l-i-l-^D , 


"12 


~ 2 2^ 2-^ 2^ ~ 2^' 


"14 = 


2 + 2^ ~ 2-^' + 2^ 2^' 


/3l2 




/9l4 = 


-J, 


7l2 




7l4 = 


-A;, 


(^12 


= -3, 


5i4 = 




012 


^-i-j-k-^D, 


014 = 


-i-k-l-\D, 


Pl2 


= l + k-l, 


PlA = 


1 + ^ - J, 


0'l2 


= —l^i — 23 ~ 2^ 2^' 


0i4 = 


~2^ 2-^ ~ 2^ ~ 2^ ■ 



Observe that when the parameters 6 and 9' in 5i are equal, then our defined 
function Si becomes the known Appel's hypergeometric function F^, whereas 
when the parameters and 0' in S2 are equal, our defined function S2 reduces 
to the known Appel's hypergeometric function H2 (see Section 3). 

Looking carefully at these one can verify without difficulty that there are 
symmetry relations among them. For example, if we make the substitution 
s t, i k and j' <-> Z in (2) we obtain (3). In a similar manner, in (4) the 
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substitution k ^ I transforms it in (5) and the substitution s t, j k, 
i ^ I yields (5) ^ (6). There are several other symmetry properties of the 
box diagram which transform one solution into another. 

Now we must combine them in such a way to have sums of linearly indepen- 
dent solutions bearing the same functional variable. This is the constructive 
prescription [1]. We then get from the above list six types of functional vari- 
ables, that is, six new such combinations or six new results for the Feynman 
integral (1), namely, 

J3 = h, J4 = I7, (14) 

J5 = /8 + /9, J6 = /l0 + /ll, (15) 

Jj^Iu + hs, J8^Iu + h5- (16) 

Note that the relevant Feynman integral is obtained via 

K{i, j, k, l;m) = J (-i, -j, -k, -I; m) 



3 Regularization and Discussion. 

Here we constrain ourselves to the special case where the integral (1) is the 
one for QED photon-photon scattering at the one-loop level, that is, we are 
interested in taking the particular values i = j = k = I = —1. However, these 
expressions become singular when we take the referred limit and/or let D — A. 
Therefore, some kind of regularization procedure is called for. 

For the first two, i.e., J3 and J4, we can adopt the standard procedure of 
dimensional regularization [6]. Introduce D = A — e and expand the whole 
expression around £ = to get 

F3(l,| + i£,l,l + i£;2 + |£|x,y),(17) 

where x = —t/s, y = 4m^/s, F3 is a hypergeometric function of two variables 
which is absolutely convergent for \x\ < 1 and \y\ < 1, and je is the Eu- 
ler's constant[10,12]. We can write a simpler expression by using a reduction 
formula [10,12,13], 



tR 
-'6 



-2 



+ log (-27rs) + 7b 
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where z = y/{y — 1) and Fi is another hypergeometric function of two vari- 
ables which is absolutely convergent in the same region of the F3 above. This 
function has a simple integral representation[12], 

F^{a,P,P';^\z„Z2) = .^V ^ r / du- ^, (19) 

r{a)r{'y — a) J [1 — uzi)p [1 — UZ2)'^ 



where the parameters must satisfy Re(Q;) > and Re(7 — a) > 0. It is straight- 
forward to evaluate this integral when the parameters take the values we have 
in hands. Substituting (19) and (18) in (17) and expanding the hypergeometric 
function in Taylor series, we get 



J3(-l,-l,-l,-l;m) 



rR 



s{s — Am?) 



7b -Flog 1 



■ dp> - + log (-27rs) 
F^{a,P,P'■n\x,z). (20) 



Note that there is a simple pole which we did not expect by naive power 
counting. We will discuss this singularity and the one that appears in the 
following solution in the next subsection. Here we introduce the parametric 
derivatives[8,14], 

= dM^) = {a\z) ma + z)- ^^;{a)] , (21) 



where the ■0— function is the logarithmic derivative of the gamma function 
[12,15]. First carry out the parametric derivatives in (19) then substitute the 
values of the parameters and integrate. For the other terms the integral results 
in, 

F,(l,l.i;2|...)^-|^log(i±f^f-|i), (22, 

I -Kst V i — Kst ris + ^st / 



where 

„ / 4m^ „ / Am? Am? , , 

Rs = ^\-^. Rst = ]Jl-- -. (23) 



See that the limit i — > is well-defined. 



We can write down immediately the result for the integral Jy by noting that 
it can be transformed into Iq if we make the changes i k, j I and s t, 
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J4(-l,-l,-l,-l;m)=/7" 



R 



t{t - Am?) V e 



+7b + log 1 - 



t 



- df}, -d^ + log (-27rt) 
F,{a,(5,(5'-^\w,w').{2A) 



where w — —s/t and w' — Am? /{Am? — t). For the region of convergence see 
figure 1. 

For the remaining solutions dimensional regularization is unsuitable to reg- 
ularize their divergences. Consider for example (4) where there is a factor 
— k + I) which is divergent in the particular limit we arc interested in, 
i.e., i = j = k = I = —1. This factor has no D-dependence and dimensional 
regularization here is useless. What we must do is to use a different procedure, 
namely, regularizing the exponent of some of the propagators[l,8,16]. 

Let us then consider the fifth solution of the Feynman integral, J5. We must 
regularize one exponent of one of the propagators, say, k = 1 — ( (we could 
also take the exponent /.). The important point is that the final result will be 
independent of this choice. The other exponents are set to minus one while 
the dimension of the space-time remains arbitrary. Doing this we have 



r(3-iD + c)r(c)r2(iD-2-C) 



gt3-D/2 ^cr(i + c)r{D - 4 - c) 



-t -Am^ 



and 



(25) 



tR 
-'9 



r(3 - \D)v\\D - 2) r(-c) 

st^-D/2 sCT{D - 4 - C) 

xH^i^, 1 + C,1 + K' HK'l + C 



-t -Am^ 



t 



(26) 



The hypergeometric function H2 is defined by the double sum [12], 

H2{a,(3,-f,6;p\x,y) 



g {a\m - n){(3\m){-f\n){S\n) x'^y'' ^^7^ 



m,n=0 



(p|m) 



mini 



The region of absolute convergence of this function H2{. . .\x,y) is bounded 
by the lines[12], 

1,1 

\y\ < T—ri' \^\ < 1' |y| < 1 > (28) 
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see figure 2. 

In proceeding our analysis of the new results, for the solution J5, let us now 
expand the H2 function in Taylor series around C = 0) keeping terms up to 
the first order in C, 



, , , , , 2'-^{-nri^V^T{?>-\D)T{\D-2) ^ 

^5(-l, -1 , -1, -1; m) = „.3-n/9T^/l 3^ [-lE 



+ logy -2V^ (lD-l)+i,{2,-\D) + ^—^-da-d^ 
-t -4m2\ 



XH2 a, 1, 1, 



t 



(29) 



where the parametric derivatives must be taken at the point a = 0; p = 1. 
The nature and meaning of these singularities will be the touched on in the 
following subsection. 

In a similar manner we regularize the sixth solution. But now we take i — 
— 1 — C- As a result we get, 



M-l, -1 , -1, -1; m) = ts^-n/.rC,D - |) ^'^^ 

+ log(^)-2^ {lD-l)+ij{3-lD) + -^-da-dp 



XH2 { a, 1, 1, l;p 




(30) 



and as we shall see later on, these two functions H2 are related to the functions 
F3 that arc divergent too. The region of convergence can be constructed as we 
did above for the H2 function (see fig. 3). 

Consider now the seventh solution of the Feynman integral, J7. Like the pre- 
ceding case, it has a simple pole in the exponents, so that it demands only 
one suitable parameter to regularize it. Looking at (8) and (9) we note that a 
good choice to introduce our regularization parameter is to take I = —1 — 
while the other exponents can be set to i = j = /c = — 1 without any problem. 
Then, 



-'12 — 



-TT 



m'^s 



---l + logs + 0(C) 



XS2 + k, 1 + C, 1, 1, 3 - iL); 1 + c, 1 - k 



Arri^ -t 
s ' 4m^ 



, (31) 
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and 



tR. 
-'13 



-l-7E-log(-m2) + 0(C) 



XcS2 -I - \C 1, 1, 1, 3 - |D + C; 1 - C, 1 + K 



4m^ -t 



• (32) 



Now expand the factors of (31), (32) and the series (13) around C = ^"^^ 
substitute the values a = —\, I3 = ^ = 5 = p = (j) = (l)'=\. Using the 
fact that dfj + dp = Q (only because these two parameters are equal) and an 
analogous relation between and 0', we get, in four dimensions. 



J7(-l,-l,-l,-l;m) 



TT 



( —s 



2 + jE + da + dp-log{ — 



xH2 a, 1, 1, l;p 



4m^ -t 



s ' 4m^ / 



(33) 



Note that the above result is finite and that there is no dependence on and 
0', SO ^2 reduces to H2. The pole cancels out and then we can take the limit 
of vanishing (. 

The next two solutions follow the same procedure, yielding 



J8(-l, -1, -1, -l;m) 



71 



2 + ^E + da + dp- log (-^ 



mH 

xi^sl a,l,l,l;p 



4m^ —s 



t ' 4m^ / ' 



(34) 



which is also finite. The region of convergence of this solution is shown in figure 
4. We do not need to calculate the parametric derivatives because Davydychev 
already did it [8]. Using the transformation formula between H2 and ^2(13], 



H2{a, P, 7, 6;p\x, y) = Ai F2 + 7, /3, 7; p, 1 + 7 - 6 

+A2 F2 I a + 5, /3, 5; p, 1 + 5 - 7 



where we define the coefficients 



r(l - a)T{6 - 7) 

r((5)r(i - q; - 7) 



y ^, ^2 = ^i(7^<^) 



(35) 



(36) 
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we can identify the parametric derivatives of H2 with the ones of F2 calculated 
by Davydychev. Care must be taken with (35) because with the particular 
parameters we have in hands the individual terms on the RHS are singular, 
but which cancels out at the end when both terms are added together. 



3.1 Discussion. 



As we have mentioned earlier, the set of new solutions we have obtained here 
contains singular solutions that deserve a closer look. Let us examine them in 
order to understand the meaning and the nature of such singularities. To begin 
with, let us give some arguments to show the correctness of our results. As we 
had conjectured in [1], solutions containing the H2 hypergeometric functions 
did in fact appear. 



Consider the first result we obtained in our previous work[l,8], i.e.. 



Ji(-l,-l,-l,-l;m) 



t 



Aw? ' Aw? J 



(37) 



The hypergeometric function F3 which appears here is related to the hyper- 
geometric function H2 via analytic continuation (see Erdelyi[13]), 



F^{a, a', P, 7|x, y) = Bi i/2 (l + « - 7, a, a', 1 + a - /3 -y^ (38) 

+B2H2 (l + /3-7,/3,a',/?';l + /3-a 



where the two coefficients are 

^ r(/3-c.)r(7) 



B2 = B^{a^(5) 



(39) 



So, with the help of equation (38) we rewrite F3 in terms of H2 without 
worrying very much about constant factors because they arrange themselves 
properly in the process. Indeed, in this case both factors on the RHS containing 
gamma functions are singular (this is a special case of analytic continuation 
known as the logarithmic case), but whose singularities cancel out at the end, 
leaving us with a finite result as it should be. Then, 



Ji ~ F3 



1,1, 1,1; 2 



t 



4m^ ' 4m^ 



) =Cii/2(-^,l,l,l;l 



Aw^ -t 



+ 
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+C2i/2(--,l,l,l;l 



s ' Am? ) ' 



(40) 



which clearly portrays the same H2 function we have in J7. Conclusion: NDIM 
provides, even if we did not know (38) a priori, the transformation Ji — > J7, 
or, in other words, the analytic continuation formula F3 — > H2. Moreover, 
as Erderlyi[13] mentioned, there is a transformation similar to (38) for the 
variable y in F3. This will give Ji ^ Jg. 

In order to verify that there are branch points in the Fcynman integral, we 
can do the following. Consider the definition of the hypergeometric function 
H2 given in (27). Substituting the values of the parameters — recall that the 
derivatives of an analytic function are also analytic having the same region of 
convergence — two of them cancel out and we get 

H2 (-i 1, 1, 1; 1 y)=t ^''"y'"^ (-1 - H/^)^ , (41) 



where we have used the identity (a| — k) = (—1)'^/ (1 — a\k). Observe that the 
series in is a hypergeometric function iFo[15] that can be summed. It results 
in the following 

^ (^1^) ^' 



with variables x and y given in either J7 or Jg. The remaining series in i/ is a 2-F1 
hypergeometric function that can be written down in terms of an elementary 
function and it is straightforward to show that it has branch points, see (46) 
below. 

The same procedure can be apphed to J3. Using (38) for the hypergeometric 
function F3, we get 



s s \ 2 



-s —Am 



t ' 



+ 



1 



+C4if2lO,l,l,-;l 



—s —Am^ 



(43) 



yielding ^ Jq. The analogous transformation for the variable y in F3 yields 
~^ J 5 a-nd so on. 



Using analogous routes we used above, it is possible to express this H2 function 
in terms of an elementary function, this time a square root. Considering its 
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definition, the cancefing of the parameters (3 and p for the specified values and 
summing the series in u we get, 



^2 (0,1,1, 1; 1 



00 



[-|/(l-,x)]^(i|i/) 



1 



(44) 



observe that the square root in the denominator is equal to Rst, see eq.(23). An 
important point to note here is that even though the results remain divergent, 
they are still connected by an analytic continuation formula. The questions 
that need to be addressed then are now: What does this mean? What is the 
nature of these singularities? 

First of all, it is known[17] that a four-point graph like the one in the photon- 
photon scattering has no leading singularities in the physical region. Such 
singularities does happen to occur in four-point functions when the two in- 
coming particles enter the same vertex and the two outgoing particles also 
leave the same vertex (see Fig. 5). Since this is not our case, we thus conclude 
that the singularities we have do not occur on the physical sheet, i.e., they are 
harmless [17]. 

Secondly, in analytically continuing a given function from a region TZi into 
another region, TZ2 it is important that no singularities be present between 
the regions, otherwise the result for the analytic continuation may not be 
unique. The non-uniqueness always manifest itself whenever the singularity is 
of the branch point type [18]. We know that for the photon-photon scattering 
process we have a branch-cut in s = 4m^ in the s-channel, so that in carrying 
out our analytic continuation from Jq — > J3, we are crossing this branch-cut, 
and then the singularities do arise. 

This naive argumentation shows us the great possibilities of NDIM. It re- 
produces three general — with no restriction in the parameters — analytic 
continuation relations between Appel's hypergeometric functions which are 
far from trivial to obtain (see [13]). It is clear too that the technique allows 
us the bonus by-product of pinpointing singularities of Feynman integrals. 

Eden [19] devised a technique to find out the singularities of integral representa- 
tions. In [17] Eden et al applied it in the general box diagram and the equation 
of Landau's surface — the surface of possible singularities of a integral rep- 
resentation — is given by a 4x4 determinant. In our case (equal mass for the 
virtual matter fields and on-shell photons) the Landau's equation [6, 17, 20] is. 




(45) 
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so that there are four possible solutions, 



Arri^s 



(46) 



s = 0, t = 0, s = 



t = 



s — Am? 



Just here it is important to observe that the two last solutions make the 
hypergeometric function H2 in (44) and in (42) singular. They are branch 
points of the Feynman integral. The first two are the so-called pseudo-threshold 
[6,17,20,21] — singularities of the Feynman integral which occur on an unphys- 
ical sheet — see also that the possible singularities of (1) are located in the 
region of convergence of the two above functions, J3 and J4. We think that 
the poles does not cancel because of this reason, the so-called pinch singulari- 
ties. We can verify, comparing the analysis contained in [6], that the last two 
solutions of the Landau's equation are in fact singularities of (1). 



4 Conclusion. 

Using the technique known as negative dimensional integration we obtained six 
new results for the massive Feynman integral for the box diagram which con- 
tributes to the photon-photon scatteringin QED. These results are expressed 
in terms of Appel's hypergeometric functions F3 and H2. All these new results 
had to be regularized. Four of them remain divergent even after the regular- 
ization procedure, but the tenacious singularities are harmless because they 
do not he in the physical region. The other two become finite after the regu- 
larization procedure and these are important in treating relativistic dynamics 
of forward scattering. We have shown that these new solutions are correctly 
related with each other by analytic continuation. To further check these re- 
sults, calculation of the forward amplitude for the photon-photon scattering 
amplitude in QED is in progress. 
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